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Abstract 

This is a pedagogical paper which explains some ideas in cosmology at a level 
accessible to undergraduate students. It does not use general relativity, but 
uses the ideas of Newtonian cosmology worked out by Milne and McCrea. 
The cosmological constant is also introduced within a Newtonian framework. 
Following standard quantization procedures the Wheeler-DeWitt equation 
in the minisuperspace approximation is derived for empty and non-empty 
universes. 



1 Introduction 



Some of the modern ideas in cosmology can be explained without the need to 
discuss general relativity [Landsberg and Evans (1977)]. The present paper 
represents an attempt to do this based on Newtonian mechanics. There is 
a need for pedagogical articles which discuss forefront ideas in research, but 
also include all the relevant derivations in one place and at an accessible 
level. This will be useful for instructors not familiar with cosmology but who 
would nevertheless like a summary in a single article that can be used in, 
say, a modern physics course. There are also some ideas presented here that 
cannot be found in the pedagogical literature. These include a discussion of 
how to incorporate the cosmological constant in Newtonian mechanics and a 
discussion of the Wheeler-DeWitt equation for flat, open and closed empty 
and non-empty universes. 



In what follows the equation of state for matter and radiation will be needed. 
In particular an expression for the rate of change of density, p, will be needed 
in terms of the density p and pressure p. (The definition x = where t is 
time, is being used.) The first law of thermodynamics is 



where U is the internal energy, W is the work and Q is the heat transfer. 
Ignoring any heat transfer and writing dW — Fdr — pdV where F is the 
force, r is the distance, p is the pressure and V is the volume, then 



Assuming that p is a relativistic energy density means that the energy is 
expressed as [Guth and Steinhardt (1989)] 



2 Equation of State 



dU + dW^ dQ 



(1) 



dU = -pdV. 



(2) 



U^pV 



(3) 



from which it follows that 



U^pV + pV^ -pV 



(4) 
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where the term on the far right hand side results from equation (Q). Writing 

V 



V (xr'^ imphes that ^ = 3-. Thus 



p = -3(p + p)- (5) 



2.1 Matter 

Writing the density of matter as 



(6) 
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it follows that 

. _ dp . r 

p=—r = -3p- 7 
ar r 

so that by comparing to equation (||), it follows that the equation of state 
for matter is 

p = 0. (8) 

This is the same as obtained from the ideal gas law for zero temperature. 
Recall that in this derivation we have not introduced any kinetic energy, so 
we are talking about zero temperature. 



2.2 Radiation 

The equation of state for radiation can be derived by considering radiation 
modes in a cavity based on analogy with a violin string [Kubo (1967)]. For 
a standing wave on a string fixed at both ends 

where L is the length of the string, A is the wavelength and n is a positive 
integer (n = 1, 2, 3 ). Radiation travels at the velocity of light, so that 

9 r 

c=fX = f— (10) 
n 

where / is the frequency. Thus substituting / = into Planck's formula 
U = huj = hf, where h is Planck's constant, gives 

nhc 1 1 /o , , 

U = — - oc V~'/'. (11) 
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Using equation the pressure becomes 



Using p = U/V, the radiation equation of state is 

P=Ip- (13) 
It is customary to combine the equations of state into the form 

P = Ip (14) 

where 7 = 1 for radiation and 7 = for matter. These equations of state are 
needed in order to discuss the radiation and matter dominated epochs which 
occur in the evolution of the Universe. 



3 Velocity and Acceleration Equations 

The equation which specifies the speed of recession is obtained by writing 
the total energy E as the sum of kinetic plus potential energy terms (and 
using M = IvrrV) [Madsen (1995), Roos (1994)] 

E = I + U = -mr — G = -mr [H pj (15) 

2 7* 2 3 

where the Hubble constant H = m is the mass of a test particle in the 
potential energy field enclosed by a gas of dust of mass M, r is the distance 
from the center of the dust to the test particle and G is Newton's constant. 
Recall that the escape velocity is just v escape = \/^^^ = \J^^pr'^, so that 
the above equation can also be written 

= vL,,e - k' (16) 

with k' = — The constant k' can either be negative, zero or positive 
corresponding to the total energy E being positive, zero or negative. For a 
particle in motion near the Earth this would correspond to the particle escap- 
ing (unbound), reaching infinity with zero speed (critical case) or returning 
(bound) to Earth because the speed r would be greater, equal to or smaller 
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than the escape speed f escape- Later this will be analagous to an open, flat 
or closed universe. Equation (|15D is re-arranged as 

H' = ^P + ^- (17) 

Defining k = — and writing the distance in terms of a scale factor a and 
a constant length s as r{t) = a{t)s, it follows that - = - and - = -, giving 
[Madsen (1995), Roos (1994)] 

a 3 

which specifies the speed of recession. The scale factor is introduced because 
in general relativity it is space itself which expands. Even though this equa- 
tion is derived for matter, it is also true for radiation. The same equation is 
obtained in general relativity [Islam (1992)]. According to Guth [Guth and 
Steinhardt (1989)], k can be rescaled so that instead of being negative, zero 
or positive it takes on the values —1,0 or +1. From a geometric, general 
relativistic point of view this corresponds to an open, fiat or closed universe. 

In elementary mechanics the speed i; of a ball dropped from a height r 
is evaluated from the conservation of energy equation as f = y/2gr, where 
g is the acceleration due to gravity. The derivation shown above is exactly 
analagous to such a calculation. Similarly the acceleration a of the ball is 
calculated as a = g from Newton's equation F = mr, where F is the force 



and the acceleration is f = The acceleration for the universe is obtained 
from Newton's equation 

_ G—^ = mr. (19) 

Again using M = ^nr^p and ^ = f gives the acceleration equation 

F r a AtxG 

— = - = - = —p. 20 

mr r a 6 

However because M = ^nr^p was used, it is clear that this acceleration 
equation holds only for matter. In our example of the falling ball instead of 
the acceleration being obtained from Newton's Law, it can also be obtained 
by taking the time derivative of the energy equation to give = ^ = = 
{^/2gr){^j2g^^) = g. Similarly, for the general case one can take the time 
derivative of equation ([T8|) (valid for matter and radiation) [Madsen (1995)] 

_i^ = 2Aa = — -(p.^). (21) 
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Upon using equation the acceleration equation is obtained as 

a 4:71 G , „ , 4:7rG ,^ , 

- = ^(P + 3p) = -^(1 + 7)P 22 

which reduces to equation (^) for the matter equation of state (7 = 0). The 
same equation is obtained in general relativity [Islam (1992), Milne (1934), 
McCrea and Milne (1934), Bondi (1961)]. 



4 Cosmological Constant 

In both Newtonian and relativistic cosmology the universe is unstable to 
gravitational collapse. Both Newton and Einstein believed that the Universe 
is static. In order to obtain this Einstein introduced a repulsive gravitational 
force, called the cosmological constant, and Newton could have done exactly 
the same thing, had he believed the universe to be finite. 

In order to obtain a possibly zero acceleration, a positive term (conven- 
tionally taken as |- ) is added to the acceleration equation (^) as 



AtcG , ^ A . . 

-(p + 3p) + - (23) 



a 3 3 

which, with the proper choice of A can give the required zero acceleration 
for a static universe. Again exactly the same equation is obtained from the 
Einstein field equations [Islam (1992)]. What has been done here is entirely 
equivalent to just adding a repulsive gravitational force in Newton's Law. 
The question now is how this repulsive force enters the energy equation (|T8|). 
Identifying the force from 



^ ^ ^repulsive A 



r a mr 3 

and using 



(24) 



3 dr 

gives the potential energy as 



_ A ^ _dU 

-f repulsive — —TflT = - \^^) 



Urepuisive- ^ ^mr^ (26) 



6 



which is just a repulsive simple harmonic oscillator. Substituting this into 
the conservation of energy equation 



E = T + U = -mr — G -—mr — -mr [H p 27 

2 r 23 2^ 3^3^^^ 



^' = (r)' = ^P-3 + -- (28) 



gives 

,a,o SttG A; A 

Let us comment on the repulsive harmonic oscillator obtained above. 
Recall one of the standard problems often assigned to students in mechanics 
courses. The problem is to imagine that a hole has been drilled from one side 
of the Earth, through the center and to the other side and to show that if a 
ball is dropped into the hole, it will execute harmonic motion. The solution 
is obtained by noting that whereas gravity is an inverse square law for point 
masses M and m separated by a distance r as given by F = G^^, yet if one 
of the masses is a continous mass distribution represented by a density then 
F — G^npmr. The force rises hnearly as the distance is increased because 
the amount of matter enclosed keeps increasing. Thus the gravitational force 
for a continuous mass distribution rises like Hooke's law and thus oscillatory 
solutions are encountered. This sheds light on our repulsive oscillator found 
above. In this case we want the gravity to be repulsive, but the cosmological 
constant acts just like the uniform matter distribution. 

Finally authors often write the cosmological constant in terms of a vac- 
uum energy density as A = SnGp^jac so that the velocity and acceleration 
equations become 



,a,9 StiG k A SnG , , k 



and 



d 4:7rG ,^ , A 471^,^ , 8nG 
a 6 3 3 3 

4.1 Einstein Static Universe 

Although we have noted that the cosmological constant provides repulsion, 
it is interesting to calculate its exact value for a static universe [Atwater 
(1974), Adler et al (1975)]. The Einstein static universe requires a — ao — 
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constant and thus d = d = 0. The case a = will be examined first. From 
equation ( ]23|) this requires that 



A = 47rG(p + 3p) = 47rG(l + 7)p. (31) 

If there is no cosmological constant (A = 0) then either p = which is an 
empty universe, or p = — |p which requires negative pressure. Both of these 
alternatives were unacceptable to Einstein and therefore he concluded that 
a cosmological constant was present, i.e. A 7^ 0. From equation (|3l|) this 
implies 

P = 4vrG(l + 7) ^''^ 



and because p is positive this requires a positive A. Substituting equation (|32D 



into equation (|28D it follows that 



A = ^l(f)^ + ^]. (33) 

Now imposing d = and assuming a matter equation of state (7 = 0) implies 
A = However the requirement that A be positive forces k = +1 giving 

A = ^ = constant. (34) 
% 

Thus the cosmological constant is not any old value but rather simply the 
inverse of the scale factor squared, where the scale factor has a fixed value 
in this static model. 



5 Conservation laws 

Just as the Maxwell equations imply the conservation of charge, so too do 
our velocity and acceleration equations imply conservation of energy. The 
energy-momentum conservation equation is derived by setting the covariant 
derivative of the energy momentum tensor equal to zero. The same result is 
achieved by taking the time derivative of equation (|^). The result is 

p + 3(p + p)- = 0. (35) 
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or 

and from equation (0), 3(p + p) = (3 + 7)p, it follows that 

|(pa3+-)=0. (37) 

Integrating this we obtain 

where A is a constant given by A = p^a'^'^^ . This shows that the density falls 
as for matter and for radiation as expected. 

Later we shall use these equations in a different form as follows. From 
equation 

p' + 3(p + p)^ = (39) 

where primes denote derivatives with respect to a, i.e. x' = dx/da. Alterna- 
tively 

^{pa^) + Spa"^ = (40) 

so that 

which is consistent with equation (E3) 



6 Quantum Cosmology 

The discussion of the Wheeler-DeWitt equation in the minisuperspace ap- 
proximation [Hartle and Hawking (1983), Kolb and Turner (1990), Atkatz 
(1994), Atkatz and Pagels (1982)] is usually restricted to closed {k = +1) and 
empty (p = 0) universes. Atkatz [Atkatz (1994)] presented a very nice dis- 
cussion for closed and empty universes. Herein we consider closed, open and 
flat and empty and non-empty universes. It is important to consider the pos- 
sible presence of matter and radiation as they might otherwise change the 
conclusions. Thus presented below is a derivation of the Wheeler-DeWitt 
equation in the minisuperspace approximation which also includes matter 
and radiation and arbitrary values of k. 
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The Lagrangian is 

with K = Note that this Lagrangian is not defined for a flat {k = 0) 
universe. Also for an open universe {k = —1) the Lagrangian is complex 
(A;3/2 = -t). 

The momentum conjugate to a is 

Substituting L and P into the Euler-Lagrange equation, P — ^ = 0, equa- 
tion (^) is recovered. (Note the calculation of ^ is simplifled by using the 
conservation equation ( |39D with equation (|14D, namely p' = —(3 + j)p/a). 
The Hamiltonian H = Pa — L is 

which has been written in terms of d to show explicitly that the Hamiltonian 
is identically zero and is not equal to the total energy as before. (Compare 
equation (p9|)). In terms of the conjugate momentum 

which, of course is also equal to zero. Making the replacement P ~'^'§^ 
and imposing = results in the Wheeler-DeWitt equation in the min- 
isuperspace approximation for arbitrary k and with matter or radiation (p 
term) included gives 

Using equation ( pHf ) the Wheeler-DeWitt equation becomes 

<- A? + wv^" - k'-r + —^^ = »■ 

This just looks like the zero energy Schrodinger equation [Kolb and Turner 
(1990)] with a potential given by 
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For the empty Universe case of no matter or radiation (A — 0) the po- 
tential U{a) is plotted in Figure 1 for the cases k = +1,-1 respectively 
corresponding to closed [Kolb and Turner (1990)] and open universes. It can 
be seen that only the closed universe case provides a potential barrier through 
which tunneling can occur. The open universe would immediately recoUapse. 
(Actually for the open universe the Lagrangian is complex, so this discussion 
may be meaningless. In either case it is clear that an open universe cannot 
arise via quantum tunneling.) This provides a clear illustration of the idea 
that only closed universes can arise through quantum tunneling [Atkatz and 
Pagels (1982)]. For a flat universe (/c = 0), the potential is not defined and 
thus this process of quantum tunneling cannot give rise to a flat universe. If 
radiation (7 = 1 and A 7^ 0) is included then only a negative constant will 
be added to the potential (because the term a^~'' will be constant for 7 = 1) 
and these conclusions about tunneling will not change. The shapes in Figure 
1 will be identical except that the whole graph will be shifted downwards by 
a constant with the inclusion of radiation. (For matter (7 = and A 0) a, 
term growing like a will be included in the potential which will only be im- 
portant for very small a and so the conclusions again will not be changed.) 
To summarize, only closed universes can arise from quantum tunneling even 
if matter or radiation is present. 

This work was supported by the Wisconsin Space Grant Consortium. 
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Figure 1 

Wheeler- De Witt potential for open and closed universe. 
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